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Abstract. We study the electron transport through the quantum dot coupled to 
the normal metal and BCS-\ike superconductor (N — QD — S) in the presence of the 
Kondo effect and Andreev scattering. The system is described by the single impurity 
Anderson model in the limit of strong on-dot interaction. We use recently proposed 
equation of motion technique for Kcldysh nonequilibrium Green's function together 
with the modified slave boson approach to study the electron transport. We derive 
formula for the current which contains various tunneling processes and apply it to study 
the transport through the system. We find that the Andreev conductance is strongly 
suppressed and there is no zero-bias (Kondo) anomaly in the differential conductance. 
We discuss effects of the particle-hole asymmetry in the electrodes as well as the 
asymmetry in the couplings. 



PACS numbers: 74.50.+r, 72.15.Qm, 73.23.Hk 



Electron transport through strongly interacting quantum dot 



2 



1. Introduction 

In recent years there has been much experimental and theoretical work on electron 
transport through nanometer-size areas (metallic or semiconducting islands) containing 
small number of electrons. These islands (sometimes called the quantum dots) are 
coupled via tunnel barriers to several external electrodes making it possible to adjust 
the current flowing through the system pQ . The devices give a new possibility of studying 
several well-known quantum phenomena in novel and highly controllable way. For 
instance, it is well known, that quantum dot behaves like magnetic impurity in a metallic 
host and in particular displays the Kondo effect 0-0 • Kondo effect is a manifestation 
of the simplest state formed by the impurity spin and conduction electron spins. This 
state gives rise to a quasiparticle peak at the Fermi energy in the dot spectral function 
jHj-[0] and zero-bias maximum in the differential conductance observed experimentally 

Another example is the Andreev scattering ^H], according to which an electron 
impinging on normal metal - superconductor interface is reflected back as a hole and 
the Cooper pair is created in superconductor. This effect has been shown to play crucial 
role in the transport properties of various hybrid mesoscopic superconducting devices 
[To] . There is a number of papers in the literature concerning the electron transport in 
various realizations of such devices. Here we are interested in study of the normal state 
quantum dot coupled to one normal and one superconducting electrode (N — QD — S). 
Such system was studied within scattering matrix technique [13 HE]. However this 
approach is valid only for noninteracting systems and cannot take into account effects of 
Coulomb interactions between electrons on the dot, which are very important in these 
small systems, as they lead e.g. to the Coulomb blockade phenomena ^H] or Kondo 
effect 0-0- Transport through noninteracting quantum dot has also been studied 
within nonequilibrium Green's function technique. The effect of multiple discrete levels 
of the dot was discussed in Refs. [2UI |2T] . the photon-assistant transport in [22| , electron 
transport in the weak magnetic field in j2Hj, temperature dependence of the resonant 
Andreev reflections in [24 and transport in three terminal system (two ferromagnetic 
and one superconducting electrodes) in |23] . 

In the presence of strong Coulomb interaction in N — QD — iV device the Kondo 
effect appears an influences the electron transport in the system. If one of the electrodes 
is superconducting both single electron and the Andreev current is affected by the 
Abrikosov - Suhl resonance. This problem has been extensively studied within various 
techniques j2H EBJ-jHSj and there is no consensus. Some authors have predicted 
suppression [2U EEJ EDI dJj of the conductance due to Andreev reflections while others 
- enhancement [23 EH]- Recently it has been shown [32] that one can obtain either 
suppression or enhancement of the conductance in dependence on the values of the 
model parameters. Recently other effects like emergence of the Kondo-like peaks in the 
local density of states (DOS) at energies equal to ±A (A is the superconducting order 
parameter) [3U1 1221 or a nove l co-tunneling process (Andreev-normal co-tunneling) [33] 
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have been revealed. This process involves Andreev tunneling from the QD — S interface 
and normal tunneling from iV — QD interface. As a result the Cooper pair directly 
participates in the formation of the spin singlet (Kondo effect) and leads to emergence 
of the additional Kondo resonances in the local DOS and enhancement of the tunneling 
current. 

The purpose of the present work is to apply the new technique to derive formula 
for the current through QD (in the limit of strong on-dot Coulomb interaction) in 
terms of various tunneling processes. We also study the interplay between the Kondo 
effect and Andreev reflections to give additional insight into the the problem of the 
suppression/enhancement of the zero-bias current- voltage anomaly. Further we discuss 
the question of participation of the superconducting electrons in creation of the Kondo 
effect. And finally we investigate the influence of the electron - hole asymmetry in the 
leads on tunneling transport as well as the asymmetry in the couplings to the leads. 

The paper is organized as follows. In section|2]we present model under consideration 
and derive formula for the current using EOM for nonequilibrium GF. In section El we 
apply the obtained formula for the current to the numerical study of the transport 
through N — QD — S system. Section H] is devoted to the interplay between the Kondo 
and Andreev scattering. In sectionElwe discuss the quantum dot asymmetrically coupled 
to the leads, while the effect of electron - hole asymmetry in the leads is investigated in 
sec. Some conclusions are given in sec. 



2. Model and formulation 



The Anderson Hamiltonian of the single impurity [33], m the Nambu representation, 
can be written in the form 



Ako- 

E^^ + E^^ (1) 
<y Ako 

where the Nambu spinors \I/Akcr and $ CT are defined as 

*«- = (T ) ; = ( t ) (2) 



c A-k-o ' v a - 



and 



Hi, = ( AAk ) (3) 

V A Ak ~ e Ak / 

denotes Hamiltonian of the normal (A^k = 0) or superconducting (Ag^ ^ 0) lead. 

h qd = ( E d + U d n_ a \ 

V -( Ed + U d n a ) ) [) 

is the dot Hamiltonian and 
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the dot-electrode hybridization. Here A = N, S denotes the normal metal (N) or 
superconducting (S) lead in the system. The parameters have the following meaning: 
c Ako- ( c Akcr) denotes creation (annihilation) operator for a conduction electron with the 
wave vector k, spin a in the lead A, A^k is the superconducting order parameter in the 
lead A (Ask = As, A^k = 0) , and VAk is the hybridization matrix element between 
conduction electron of energy eAk in the lead A and localized electron on the dot with 
the energy E&. <i+ (d a ) is the creation (annihilation) operator for an electron on the dot 
and Ud is the on-dot Coulomb repulsion. 

To derive the formula for the average current in the system we start from the 
time derivative of the charge (for convenience we perform calculations in the normal 
electrode) |3o] : 

J=-ej t (N N ) = j([N N ,H]_) (6) 

, where Nn = J2ka c Nka c ^ is the total electron number operator in the lead N and e 
is the elementary charge. The above formula can be written in terms of the Green's 
functions (GF): 

^-TEr^M M M]n (7) 

fl , J —oo Z7T 

, where G^f kda (u) is the Fourier transform of the Keldysh matrix Green's function 
G^f kd(T {t) = i($+(0) (g) ^Nka(t))- Now we have to calculate the Green's function 
^wkdffH- O ne can do this in the usual way, i. e.. using Keldysh equation [^EHI 

G<{lo) = (1 + G r {Lo)TT{u))G<{Lo) 

x(l + Y> a (u)G a (uj)) + G r (cj)S < (w)G a (cj) (8) 

(superscripts r, a are for retarded and advanced GF respectively) and make more 
or less justified approximations for the 'lesser' self-energy S < (c<j). Usually one uses 
approximation due to Ng j^Zj which states that full 'lesser' self-energy is proportional 
to the noninteracting one (S < (c<j) oc Sq (lv)). This approximation is widely used in 
the literature 12E]- However we wish to use another approach based on recently 
proposed equation of motion technique (EOM) for nonequilibrium Green's functions 
[55] . The usual equation of motion derived from Heisenberg equation yields undefined 
singularities, which depend on the initial conditions. The advantage of this new 
technique, based on Schwinger - Keldysh perturbation formalism, is that it explicitly 
determine these singular terms. Moreover, together with EOM for retarded (advanced) 
Green's functions it allows to treat the problem in very consistent way making similar 
approximations in the decoupling procedure for all types of the Green's functions. Such 
approach was recently proposed to calculate the charge on the quantum dot upon 
nonequilibrium conditions [39J. 

According to the Ref. [HB| the equation for the 'lesser' Green's functions reads: 

{(A\B))< = g<(u)([A,B\ ± ) + 

ff »(([A^]-|S»<+ 5 <(o;)(([Aiy J ]_|S))- (9) 
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where g ri<) (uj) is the free electron GF and Hj denotes interacting part of the 
Hamiltonian. In general this equation allows to calculate the the GF on the left hand 
side, however in practice we have to approximate the higher order GF-s appearing on 
the right hand side (as in EOM for equilibrium GF-s). But performing analogous 
approximations in the decoupling procedure for both retarded and 'lesser' GF we make 
theory consistent. 

Applying above equation for the Green's function occurring in the formula for the 
current (J2J) G^ kad (uj), we get following expression: 

^ k M#k(($ CT |$+>>° (10) 

where as an interacting part of Hamiltonian Hj we have taken the third term in the 
Eq.fJIJ. 9nu ( u ) i s retarded ('lesser') free-electron matrix GF of the normal-state 
electrode. 

^ k H=(^f^ i ) (ii) 

\ uj+e Nk +iO ) 

9Nk{UJ) ~ { 2mf(u + eV)5(u + e Nk ) ) {U} 

where f(uS) is the Fermi distribution function and eV = Hn — Us corresponds to the 
applied voltage between normal state electrode with the chemical potential fi^ and 
superconducting one with fis- In the following we fix the chemical potential of the SC 
electrode ([is = 0) and use eV as a measure of the bias voltage. 

Expression (fTTH) is general for Anderson model and doesn't depend explicitly on 
the form of the Hamiltonian describing quantum dot (Hqd). The dependence enters 
only through the Green's function G^ a \uj) = ((^al^t))^^ ■ m ^ ne following we wish 
to study quantum dot in the limit of strong on-dot Coulomb repulsion (Ud — > oo). In 
this limit double occupancy of the dot is forbidden and it is convenient to work in the 
slave boson representation in which the real electron operator is replaced by product 
of fermion and boson ones (d a =^ b + f a ) jlUJE]. Additionally the fact that there is no 
double occupancy on the dot should be taken into account in some way. Usually such 
constraint is added to the Hamiltonian by the Lagrange multiplier. There is a number 
of variants of this approach in the literature and here we shall work within La Guillou 
- Ragoucy scheme |12l ESI • l n this approach the constraint of no double occupancy is 
enforced through modification of the commutation relations of both fermion and boson 
operators in comparison to the standard ones. This approach was successfully used in 
the study of the charge on the quantum dot [39J. 

Hamiltonian of the system in the limit Ud —>■ oo in the slave boson representation 
is given in the form but now with 
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and 

jjQD _ I - ) (14) 




Having introduced slave boson representation, we can begin calculations of the 
advanced and 'lesser' on-dot GEs appearing in Eq. JTUJ). To do this we apply 
formula together with usual prescription for the advanced (retarded) GF. One 
can investigate that on the higher-order GFs appeared in the this process have similar 
form in both cases: 'lesser' and advanced. So idea is to make the same approximations 
in the procedure of decoupling of the higher order GFs. Explicitly, we have performed 
decoupling 

((c^cy^AlB))^ « cW<W™ Ak ((A|B» w (15) 

and neglected the other GFs. In above formula n\k is the concentration of the electrons 
in the lead A in state k and for superconducting electrode is given by 



(16) 



with quasiparticle spectrum F| k = e| k + A|, while for the normal lead this relation 
reduces to 

nNk = /(ejvk - eV) (17) 

We want to stress here, that we haven't used factorization like 

{{4^4wA\B)). ~ 5 xx ,5^{c^ct_J){{A\B)) u (18) 

The reason comes from the requirement of the hermicity relation between retarded and 
advanced Green's function, i. e.. G r (u>) = [G a (u;)]t. If we calculate retarded GF within 
EOM and perform the same decoupling as in advanced Green's function keeping also 
( c Ako- c A-ko-) terms, we get expressions for the GFs which violates the hermicity relation. 
The only way to fulfill that at this level is to make approximations due to (|15|) and 
neglect the remaining higher order GFs. 

The resulting advanced on-dot GF G°j CT (u;) can be written in the form of the Dyson 
equation: 

G3» = g a M + ^(u,)£S»G^(u;) (19) 

where g^i*^) non-perturbed dot's advanced Green's function: 

/ l-n-g q \ 

= "~ E *- i0 (20) 

V u+E d -iO J 

and self energy T,^ a (u) which can be written as sum of the noninteracting £° a (co>) and 
interacting Y> d a (uj) part 

E5(«) = + Ej» = £ + (21) 

Ak 
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+ 



(1— n_ CT ) I cj-Es k -jO a;+f? Sk -iO 



\V S 



n CT ) V 



-iO 



^--^Sk-jO 



(1 

|Vsk| 2 



"fik^Sk "Sk^Sk 

u+Esk—iO u—Esk—iO 



(1— n CT ) \ a;— i? 



Sk 



_ J ^Sk I 

-iO ^ cu+E sk -iO J ) 



(22) 



where we have introduced BCS factors w| k 
normal state corresponding expression is: 



ifua 



Id 



-gH. For the 



IVjvkj 



1 



-iOa 
J Afk 



(1— n_ CT ) aj— £]vk— * 









\v Nk \ 2 



, (23) 

(1— n CT ) (jJ+tNk— iO 

It turns out that within the present approach the interacting part of the self energy 
is simply related to the noninteracting one. Moreover the same relation also holds for 
retarded as well as 'lesser' GFs. This is a result of the consistency of the decoupling 
procedure and requirement of the hermicity relation between retarded and advanced 
GF. In general this relation can be written as: 

S AkM = ^Akr 3 S^ k (^)r 3 

where r 3 = ^ * _^ 

of the wave vector k in the lead A given by (jlfij) and (|17jl. 

It is possible to write Eq. for the 'lesser' GF in the form of the Keldysh equation 
(JHJ) with G a da {uj) given by Eq. (fTH|) and G r = Free electron dot's 'lesser' Green's 

function is given in the form: 

V 2tu(1 - n a )f(u)5j ) 

where = 5(a> ± i^) . As we have mentioned, the 'lesser' self energy has the same form 
as advanced one: 



(24) 



is the Pauli matrix and nx^ is the concentration of the electrons 



OJ 



A 



(26) 



where noninteracting part due to SC lead is: 



:° s <(uj) = 2nif(u) 



V (ftty^k^k (si 

and for the normal lead we have 



\V S k\ : 
(l-n CT ) 



(27) 



IViVkl 2 



where <5t = ± -Esk) an d 



/(w - eV)5 




N 







(28) 



(i 



5 — uyju _i_ j^sk; ciiiu. — 5(cj ± e^vk). And again, the interacting part of the 
lesser' self energy is related to the noninteracting one simply through Eq. (J24jh 
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Now we are ready to write the expression for the current (J7|) in terms of known 
GFs. First, let's rewrite the Keldysh equation for the element 11 of the dot's GF in 
the form: 

r~i< r ir v< r^ a _l r ,r v< r^ a i 

^11 — Lr ll 2j ll lj 'll + ( - t 11 2j 12 Lt 21 + 

^12^21^11 + ^12^22^21 (29) 

Note that we don't have term proportional to g£ as it vanishes in our case 

To calculate the GF given by Eq. fjlOj) . entering to the expression for the current 
(JJJ), we need yet element 11 of the advanced GF, more precisely imaginary part of that. 
Note that g^ b f k is purely imaginary and we need real part of g^ifid- We can write 
down equation for the imaginary part of the element 11 of G a d in the similar fashion as 
Eq.pty, i.e. 

ImG a n = G r u Imi: a u G a u + G r u Imi: a 12 G a 21 + 

G\ 2 ImJ: a 21 G a u + G\ 2 lmL a 22 G a 21 (30) 

Substituting now the Eqs. (f27I|) and (|3T)J) into (fTUj) we get expression for the Gjvk.d, 
which determines the current flj|). Finally the current (JZJ) can be written as 

J = Ju + J22 + J12 + J a (31) 
The first term represents conventional tunneling and is given in the form 



2e f 00 duj 

r N p»[f(u-eV)-f(uj)] (32) 



where is the elastic rate defined as = 27rV^p JV (0) and p N (0) is the bare normal state 
density of states at the Fermi energy. The second term describes the 'branch crossing' 
process (process with crossing through the Fermi surface) in the language of the BTK 
theory (Blonder - Tinkham - Klapwijk) [13]: electron from the normal lead is converted 
into the hole like in the SC lead. 

2e f°° du 3 - 



22 — Z" / — ^'1^22 ^12 

n J-00 In 

r Nf £\f(u>-eV)-f(u>)] (33) 

The next term corresponds to the process in which electron tunnels into SC picking up 
the quasiparticle and creating a Cooper pair. 

r N p»[f(oj-eV)-f(u)] (34) 

The last term in (pTTj) represents Andreev tunneling in which electron from the normal 
lead is reflected back as a hole and Cooper pair is created in superconducting electrode. 
= _2e ^ /roE „ |G 



A = — 7~ / — Jmz J22 |Uri2| 
h J-00 Ztt 

JV 



T N p» x [f{u-eV)-f{LO + eV)\ (35) 
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As one can see, at energies \eV\ < As and zero temperature, the only process which 
contribute to the total current, is the Andreev tunneling. The remaining ones represent 
'single particle' processes which are suppressed at \eV\ < As due to the lack of the 
states in superconductor. Of course for energies \eV\ > A s all these processes give rise 
to the current, even Andreev does, however is strongly suppressed, but still finite. It 
is also worthwhile to note that all these processes (except J n ) proceed through virtual 
states on the dot. 

3. Density of states 

In the following sections we will present numerical results of electron tunneling in the 
iV — QD — S system and show how different terms of Eq. (pTTj) contribute to the total 
current and differential conductance. But firstly we want to turn our attention to the 
density of states as it gives a lot of information about system. 

The most pronounced fingerprint of the Kondo effect in the N — QDN system is the 
Abrikosov - Suhl or Kondo resonance at the Fermi level and its temperature dependence. 
Kondo resonance appears as temperature is lower than parameter dependent Kondo 
temperature Tk- In the original Kondo effect there is odd number of electrons on the 
dot, so the total spin is half - integer. In this case electrons from the leads with energy 
close to Fermi level screen the spin on the dot producing resonance at the Fermi energy. 
If electrodes are made superconducting situation is more complicated, as there enters 
another energy scale - superconducting transition temperature T c (or equivalently SC 
order parameter A). And the Kondo effect takes place provided Tk > T c , otherwise is 
absent due to lack of the low energy states in the leads to screen the spin on the dot. 
Naturally there raises a question what will happen if one of the lead is superconducting 
and another in the normal state. This was investigated in [2E]-[2B], CHI and it has 
been found that Kondo effect survives in the presence of superconductivity in one of the 
electrodes, even Tk < T c . The reason for this is simply: the spin of the dot is screened 
by electrons in the normal lead. We will show that this is a really the case and this is 
seen in the density of states of quantum dot. 

In the FigG]we show the density of states of the quantum dot for various positions 
of the dot energy level It is clearly seen that Kondo effect, which manifests itself 
in the resonance on the Fermi level, survives the presence of superconductivity in one 
electrode. The additional structure at u = ±A coming from the SC lead is also visible. 
This is simply reflection of the SC gap. At this point it worthwhile to note, that if 
\Ed\ > A there is bound-like (Andreev) state within the SC gap, position of which 
depends on Ed- In the S — QD — S system this is a true bound state. However in the 
present case, due to finite DOS in the normal lead, this state acquires a finite width 
(resonance state). 

It is very interesting to see how the DOS will be look like in the nonequilibrium 
situation (eV = [An — l^s =¥" 0)- Let us recall that the Kondo resonance is located at 
the Fermi level of the lead. In the iV — QD — N system when eV ^ there emerge two 
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-0.3 -0.2 -0.1 0.1 0.2 



Figure 1. The density of states of the quantum dot for various values of the dot 
energy level Ed- Other parameters are following: Tjv = Ts = 0.02, A = 0.1, eV = 0, 
T = 10~ 5 in the units of the bandwidth W. 



resonances at Fermi levels of the left and right lead respectively. In our case there is 
a gap in the SC lead DOS, and if our simple picture that Kondo effect is only due to 
normal lead, we expect only one resonance pinned to the normal metal electrode Fermi 
level (hn)- As we can learn from the FigElthis is really the case - the Kondo resonance 
follows the Fermi level of the normal lead. 




-0.3 -0.2 -0.1 0.1 0.2 



Figure 2. Equilibrium (eV = 0) and nonequilibrium (eV = ±0.01) density of states of 
the quantum dot. Other parameters have following values: Tjv = Tg = 0.02, A = 0.1, 
E d = -0.08, T = 10~ 5 in units of the bandwidth W. 

In the Fig|H]the DOS is plotted for a few values of the SC order parameter A. As 
we can see the Kondo resonance is strongly suppressed in comparison to the N — QD — N 
DOS (solid line). The reason for that is that due to lack of the low laying states in SC, 
the spin on the dot is weakly screened. Similar conclusions have been reached by A. A. 
Clerk and coworkers jSD] within NCA approach. In the inset the density of states at 
the Fermi level p(sf) is plotted as a function of the SC order parameter A. 
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Figure 3. Equilibrium density of states of the quantum dot for various values of the 
order parameter (A). Ed = —0.2 and other parameters have values as in the Fig.©. 



Finally we want to discuss the temperature dependence of the dot density of states. 
In the Fig 0] we show the DOS for number of temperatures. As we expected the Kondo 
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Figure 4. Temperature dependence of the DOS for T N = 0.05, T s = 0.02, A = 0.05, 
E d = -0.2, eV = 0. 

resonance disappears as temperature is raised. However it is important to stress out 
that resonances and dips to ~ A are also temperature dependent. It is clearly seen in 
the FigEl where height of the peaks (a) and dips (b) are shown. The spectral weight of 
the Kondo resonance is also shown for comparison. It is worthwhile to note that below 
Tk the height of both peaks and dips are constant. As soon as temperature exceeds Tk 
height of these resonances starts to raise as well as dips does. This effect can be though 
of as a transfer of the spectral weight between Kondo resonance and Andreev states. 
Let us remind that Andreev processes still take place at T > Tk but less than T c . We 
want to notice that our results are in contradiction to what has been found in Ref.[5U] 
within NCA. The authors of jHD] have shown that the resonances at ±A do not appear 
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Figure 5. Temperature dependence of height of the peaks (a) and dips (6) near ±A 
(see Fig.(j2l). For comparison the temperature dependence of the Kondo resonance is 
also shown. The parameters are the same as in the Fig(j2Jl. 



at higher temperatures. This has a important consequences. This means that, even for 
T c > Tk, superconducting electrons do participate in the Kondo effect. So our simple 
picture breaks down. 

As we mentioned, they applied NCA technique to calculate the dot DOS, which 
is known to give correct results in the normal state for a wide range of temperatures. 
It is also known, that EOM gives quantitatively incorrect results at low temperatures 
(T < Tk)- But here this effect certainly take place at T > Tk in the range of validity 
of EOM. So the question of which picture is indeed realized in the real system seems 
to remain open. 

4. Andreev reflections and the Kondo effect 

We have shown, that Kondo peak in the density of states survives the presence of 
the superconductivity, however should we expect that peak in the current - voltage 
characteristic (differential conductance - G(eV) = dJ/d(eV)) ? If we consider tunneling 
processes, described by Eqs. (|3T |) - (j35j) . we might expect Kondo peak only in the 
tunneling mediated by the Andreev reflections (J35|) . The amplitudes of the other 
processes is equal to zero (at T = 0) for energies less than SC gap. 
Let's rewrite Eq. into the form: 

J a = ~ ^n s {u>)[f{u - eV) - f(u + eV)\ (36) 

We have introduced 'transmittance' T$ s (lj), associated with the Andreev tunneling, 
defined as: 

z#sM = |r^f>) P £H|G 21 M| 2 (37) 

In fact, at zero temperature and at energies less than superconducting gap T^ s {uj) can 
be regarded as a total transmittance, because Andreev tunneling is only process allowed 
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in these circumstances. T$ s (uj) for different values of the eV is plotted in the Fig. El 
The broad resonances at uj « ±0.06 are reflections of the dot energy level = 0.08 

4e-4 



3e-4 



In 2e-4 
<-z. 
t— 

1e-4 


-0.08 -0.04 0.04 0.08 

m 

Figure 6. T^ s (uj) for different values of the bias voltage eV = (solid line), —0.03 
(dashed) and 0.03 (dotted line). E d = -0.08, T N = T S = 0.01 and A = 0.1. Inset: 
large scale view of the equilibrium T^ s (ui). 

for electrons and holes shifted from its original position due to renormalization 
caused by the strong Coulomb interaction. But more important point is that there is 
no Kondo peak in equilibrium (eV = 0) transmittance. This is in agreement with Refs. 
[213 EDI- This is because the imaginary part of the anomalous Green's function G 12 (o;) 
behaves like \u\ for uj <C A while its real part is proportional to ujln{uj) and both vanish 
for uj — 0. And this is sufficient to suppress the Kondo effect. However as soon as 
we go away from the eV = 0, we can observe the Kondo peaks at energies uj = ±eV 
with approximately equal spectral weight. However there is strong asymmetry between 
negative (dashed line) and positive (dotted line) voltages. While in former case we have 
very well resolved resonances, in the later these resonances are strongly suppressed. 
This asymmetry is strictly related to the density of states (see Fig. (J2J)), where we also 
observe such asymmetry, which is associated with different conditions for the Kondo 
effect in both cases (note quantity E d — eV). The fact, that we observe the Kondo 
peak for both electrons (eV) and for holes (— eV) is in contradiction to what has been 
observed in Ref. [2H], where only small kink has emerged for uj = —eV. This is certainly 
due to different approximation scheme used in calculations. 

Since equilibrium transmittance T^ s {uj) doesn't show the Kondo peak, we cannot 
expect it in the differential conductance Ga{^Vsd) — oIJa/ d{eVsD) with J a defined by 
(JHnj), since Ga is proportional to the equilibrium T$ s . Indeed this is what we observe 
in the (Fig. Ej). We see, that Ga{gVsd) is very sensitive to the position of the dot 
energy level. The larger (negative) the smaller conductance. It can be understood 
as follows. The probability of the Andreev reflections depends on the density of states 
(for uj < A) of the normal electrode as well as dot itself. The later one is strongly 
-E^-dependent (see Fig. UJ). For Ed ^ there are no states on the dot participated in 
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Figure 7. The Andreev differential conductance Ga(bVsd) = AJa/ d(eVsD) for 
different values of the dot energy level Ed = —0.08 (solid line), —0.14 (dashed) and 0.2 
(dotted line). IV = T s = 0.01 and A = 0.1. 

tunneling between normal electrode and the superconductor. Note that in fact Andreev 
reflections take place between SC electrode and the dot. 

The lack of the peak in the differential conductance confirms that the Kondo effect 
is suppressed in the N — QD — S system with strong on-dot Coulomb repulsion. This 
result is in full agreement with those of Refs. HUB E2] 

5. Asymmetric coupling 

The quantum dot asymmetrically coupled to the normal state electrodes shows 
anomalous Kondo effect [13], which has also been observed experimentally ^21 EI]- 
This anomaly features in the non-zero position of the Kondo resonance in the differential 
conductance. In other words, if we increase one of the couplings to the leads (Tl(R)), 
the zero-bias anomaly moves to non-zero voltages [45]. In the N — QD — S system, there 
is no Kondo resonance in the differential conductance. Similarly there is no it in the 
equilibrium transmittance. On the other hand the Kondo peak emerges when system 
is in nonequilibrium. So in fact we could expect the non-zero bias Kondo peak in the 
differential conductance of the dot asymmetrically coupled to the leads. 

We have calculated Andreev transmittance and differential conductance for a 
number of couplings to the leads. The example is shown in the Fig. |H1 Unfortunately we 
haven't observed the Kondo peak at non-zero voltages regardless how big the asymmetry 
(Tn/Ts) was. The reason for this might be that for N — QD — N system the shift of 
the Kondo peak to non-zero values is very small [13], and in the present case the Kondo 
resonance cannot develop because there are too few states in the transmittance spectra 
around the Fermi energy. 

However the asymmetry in the couplings lead to another interesting behavior of the 
Andreev transmittance as well as differential conductance. Namely it turns out that 
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Figure 8. The Andreev transmittance T$ s for different values of the couplings to the 
leads T N = T s = 0.01 (solid line), T N = 0.015, T s = 0.01 (dashed) and T N = 0.01, 
T s = 0.015 (dotted line). E d = -0.08 and A = 0.1. 

more influences the T$ s (Ga) around energies close to value of SC gap than Ts does 
(see dashed line in the Fig. [SJ. On the other hand the transmittance (conductance) 
around l-E^I < A is more affected by Ts (compare dotted line in the Fig. |BJ). We can also 
note that the positions of the broad resonances, corresponding to the dot energy level 
Ed, depend on asymmetry This is due to the real part of self-energy. More important is 
that Tn and Ts shift the positions of the resonances in opposite directions: T^ towards 
Fermi energy fi^ = /4s = while to higher energies. 

In the present approach the Andreev transmittance (and differential conductance) 
vanishes at zero energy. However it shows interesting properties at the other 
characteristic energies of the system, like Ed or A. The Andreev transmittance T$ s 
at these energies is shown in the Fig. El as a function of the asymmetry in the coupling 
rjv/Ts. We see that the tunneling due to the Andreev processes at energies ±E d is 
likely to take a place when the couplings are more or less symmetric, i. e. T^ = Ts- 
More surprising result is that the largest probability of these processes at energies ±A 
is for large asymmetry T^/Ts ~ 10. As we already mentioned such behavior can be 
explained by renormalization of the dot energy level due to the real part of self-energy, 
which depends on T^ and in rather a complicated way . 

6. Particle-hole asymmetry 

Until now we have presented results for the special case, namely the electron-hole (e — h) 
symmetry in the leads. It is well known [2^] that the particle-hole asymmetry in the 
normal metal/superconductor tunnel junctions and metallic contacts suppresses the 
Andreev reflections due to the fact that the reflection and transmission probabilities are 
different for incident electrons and holes. In the quantum dot coupled to the normal 
and superconducting electrode we could also expect that this asymmetry will play a 
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Figure 9. T^ s at energies: —Ed (solid line), Ed (dashed), —A (dotted) and A 
(dot-dashed line) as function of the Tn/Ts for Ed — —0.08, A = 0.1 and fixed 
T N +T S = 0.02. 

role. On the other hand, this asymmetry is already present in the strongly interacting 
quantum dot (U = oo) in the Kondo regime (Ed < 0), as one studied here. However 
this asymmetry in the leads can further modify the Andreev tunneling. 

Let's start from the effect of the e — h asymmetry on the density of states. As we can 
read from the Fig. the asymmetry plays rather a minor role. The most pronounced 
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Figure 10. Density of states in various realizations of the electron-hole asymmetry 
indicated in the figure. Parameters have following values: Ed = —0.08, A = 0.01, 

rjv = r s = o.oi. 

difference is when the concentration in both N and S electrodes are changed. If change 
the concentration in one of the electrode, the effect is smaller. Moreover, the density of 
states almost does not depend on in which electrode concentration is changed. In other 
words, it seems to be sensitive to average concentration in both electrodes. 

In the Fig. ^2 we have shown the height of the Kondo resonance when electron 
concentration is varied. As one can read from the Fig. El the spectral weight of the 
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Figure 11. The height of the Kondo peak as a function of the concentrations of 
the electrons tin ( n s) hi the normal (superconducting) lead. Curve N,S - electron 
concentration in both electrodes is changed, N - in the normal electrode only and S - 
in superconducting. The parameters are the same as in Fig. 

Kondo peak strongly depends on the electron concentration in both electrodes. Moreover 
there is a strong asymmetry with respect to the n — 1 point, i. e. the peak is higher 
when the concentration of electrons is higher. It is rather expected result, as in the 
original Kondo effect the resonance at zero energy emerges due to the screening of the 
conduction electrons. So one could expect that it should depend on their concentration, 
as it does. Similar effect we observe when the electron concentration is changed in one 
lead only. It almost does not depend on in which lead n is changed. However for large 
e — h asymmetry, concentration of electrons in the normal lead seems to play a more 
important role. 

Now let's turn to the Andreev reflections and their modifications due to the e — h 
asymmetry. The Andreev transmittance T^ s (u) (Eq. (|37)l). shown in the Fig. [T1Z1 is 
also affected by the concentration of the electrons in the leads. However quantitative 
behavior of T^ s {u) seems to not depend on the e — h asymmetry. The most pronounced 
qualitative differences occur for the energies u = \Ed\. Decreasing the number of 
electrons in the normal lead, T^ s (uj) also decreases around these energies. On the 
other hand, the effect is just opposite if we decrease the electron concentration in the 
SC lead. One can rather easily explain the dependence of the T^ s (cu) on the number of 
normal electrons, taking note of the fact that probability of Andreev reflections is larger 
when number of electrons in N lead is large and the number of holes small. However 
T^ s (uj) also depends on concentration of electrons and holes in the SC lead. This can 
be understood as follow: the number of electron-like (hole-like) quasiparticles in SC is 
proportional to the concentration of the electrons (holes) in this lead in the normal state. 
In the Andreev process, if two electrons enter the SC, the electron-like quasiparticle is 
created. This means that probability of Andreev reflections depends on the number of 
electron-like and hole-like quasiparticles in SC lead. So increasing the number of holes 
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Figure 12. Andreev transmittan.ee for various electron concentrations in the leads. 
Model parameters are the same as in Fig. 



in SC, the probability of Andreev reflections of the impinging electrons is larger. This 
is exactly what we can read from the Fig. 

The modifications of the T^ s (u) due to e — h asymmetry are not so large as the 
modifications due to the asymmetry in the couplings (see Fig. |HJ), nevertheless e — h 
asymmetry influences the Andreev tunneling. 

7. Conclusions 

In conclusion we have studied a strongly interacting quantum dot connected to the 
normal and superconducting leads. Using the equation of motion technique for the 
nonequilibrium Green's functions, we derived the formula for the current in terms of 
various tunneling processes. This technique allowed us to calculate at once all the 
Green's functions emerging in the problem and perform consistent decoupling procedure 
for the higher order Green's functions. 

We discussed the problem of the interplay between Kondo effect and Andreev 
reflections. While the Kondo resonance is present in density of states, there is no 
zero bias anomaly in the differential conductance. As a matter of fact, the Andreev 
conductance is strongly suppressed for zero-bias voltages. We also further raised a 
question regarding the participation of the superconducting electrons in the Kondo 
effect. The obtained results seem to support the scenario in which they do not participate 
in the Kondo effect. 

Finally, we discussed the problem of asymmetry in the couplings to the leads and 
found the large modifications of the Andreev conductance due to this effect, mainly for 
energies around dot level and superconducting gap. We also studied the properties of 
the system when the concentration of electrons in the leads can be changed. However, 
the modifications of the Andreev tunneling due to this effect are much smaller and 
quantitative only. 
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